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We report on Landauer-Biittiker studies of anomalous Hall transport in a two-dimensional electron 
gas with Rashba spin-orbit coupling and a magnetization provided by localized magnetic moments. 
Our system is described by a discretized tight-binding model in a four-terminal geometry. We 
consider both the case of magnetically disordered systems as well as ballistic transport in disorder- 
free systems with spatially homogeneous magnetization. In the latter case we investigate both 
out-of-plane and in-plane magnetizations. We numerically establish a close connection between 
singularities in the density of states and peaks in the Hall conductance close to the lower band 
edge. Consistent with previous theoretical studies based on diagrammatic perturbation expansions, 
these peaks occur at Fermi energies where only the lower dispersion branch is occupied. Moreover, 
for large magnetization the Hall conductance is, along with the density of states, suppressed. This 
numerical finding can be understood from analytical properties of the underlying model in the limit 
of an infinite system. 
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I. INTRODUCTION 

The anomalous Hall effect (AHE) is the subject of 
a long-standing and partially still ongoing theoretical 
debateA^i^^. It amounts in a Hall conductivity 
which is not due to an external magnetic field but the 
result of the magnetization of a solid. A large portion 
of the renewed interest in this phenomenon is generated 
by research on ferromagnetic semiconductors 8 -^. In gen- 
eral, it is common to distinguish between two types of 
mechanisms for the AHE, both relying on spin-orbit in- 
teraction: The extrinsic mechanism requires the presence 
of impurities or other imperfections and is based on con- 
tributions to spin-orbit coupling from such scattering po- 
tentials. This spin dependence of the effective scattering 
potential gives rise to the skew-scattering 2 - and the side- 
jump^ contribution to the anomalous Hall conductivity. 
The intrinsic mechanism is independent of scattering cen- 
ters and is a result of the spin-orbit-coupled electronic 
band structure, where the spin-orbit interaction stems 
from the ordered crystal lattice itself. 

Among many different systems, the case of a 
two-dimensional semiconductor electron gas with 
an intrinsic effective spin-orbit coupling of the 
Rashba typei^ has attracted considerable theoreti- 
cal intereS^l ^^^ ^ ^ 8 . ^ 20 ' 21 ' 22 ' 23 and wag alsQ 

studied experimentally in an n-doped II- VI semicon- 
ductor hcterostructure containing manganese ions2^. 
The theoretical investigations have considered both the 
intrinsic effec t 17 ' 18 as well as combinations of intrinsic 
and extrinsic mechanismsi^^^^^^ 2 ^' 2 ^. In 
the present work we shall concentrate on the purely 
intrinsic AHE. 

An important tool for the theoretical description of 
transport in such mesoscopic systems is given by the 
Landauer-Biittiker formalis m 25 ' 26 . In this paper we re- 
port on numerical studies within this approach on Hall 
transport in a two-dimensional electron gas with Rashba 
spin-orbit interaction and magnetic impurities. Such an 
investigation of the AHE seems to be missing in the previ- 



ous literature although several studies of this kind on the 
related spin Hall effect are available^L 2 ^^ 3 ^^^. 

This paper is organized as follows. In section |H] we de- 
scribe important properties of our underlying model and 
outline the Landauer-Biittiker approach to transport in 
such systems. More specific information on Green's func- 
tion in semi-infinite leads used in our study can be found 
in appendix [A] In section IIIII we present our numerical 
results covering both ballistic transport and transport in 
magnetically disordered systems. We close with conclu- 
sions and an outlook in section ITVl 



II. MODEL AND APPROACH 

We consider a two-dimensional gas of non-interacting 
electrons with spin-orbit interaction of the Rashba 
typei^. Additionally, the electron spin is coupled to mag- 
netic impurities. 



A. Continuum model 

The generic single-particle Hamiltonian for the contin- 
uum system reads 
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2m* 



-r {p x cr y - PyCT x 
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(1) 



Here p is the electron momentum, m* its effective band 
mass, and the Pauli matrices a describe the electron spin. 
The strength of the spin-orbit interaction is described 
by the Rashba parameter a, and A is the effective Zee- 
man splitting due to the coupling of the electron spin 
to magnetic impurities. In general, this quantity will be 
position-dependent, A = A(r). However, it is also in- 
structive to consider the case of spatially constant mag- 
netization corresponding to an uniform impurity polar- 
ization. In the following we will consider both the case of 
a homogeneous magnetization perpendicular to the plane 



of the electron g ctS, ciS well systems with in-plane magne- 
tization. We note that an in-plane magnetization can also 
be interpreted as a genuine magnetic field which couples 
in a strictly two-dimensional situation only to the spin of 
the electron but not to its orbital degrees of freedom. For 
a homogeneous magnetization the eigenstates are given 
by plane waves and the energy of a given wave vector 
k = p/h reads 



f^k 2 
2m 

± 



e{k) 



A x ) 2 + (ak x + A y ) 2 + A 2 . (2) 



Let us first discuss the case of a purely perpendicu- 
lar magnetization, A = (0,0, A). Here one finds that, 
provided that the energy scale of the Rashba coupling 
en := m*a /h 2 is larger than the Zecman coupling, 



e R > |A| 



the lower dispersion branch in Eq. ^ has a minimum at 
finite k = kyr^in, , 




FIG. 1: Energy dispersion e±(k) for A = 0.2er (solid lines) 
and A = 1.5er (interrupted lines). Left panel: perpendicular 
magnetization A = (0,0, A). Right panel: in-plane magneti- 
(3) zation A = (A x , A y , 0) with the wave vector k being orthog- 



onal to A, A • k = 0. The energies are given in units of the 
Rashba energy er while the wave vector is measured in units 
of the inverse "Rashba wave length" kn = m*a/h 2 . 
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This dispersion minimum at finite wave vector leads to a 
van Hove singularity in the electronic density of states at 
the bottom of the lower branch, e — > e m i„ + 0. Explicitly, 
the density of states is given by 







D(s) = 



7vh 2 



2(e- 



2-rr/i 2 



£min < £ < -|A| 

-|A|<e<|A| 
|A| <e 



(6) 

This quantity has obviously a square-root singularity at 
£ — * £min + 0. As we shall see below, such singulari- 
ties are intimately linked to the observation of anoma- 
lous Hall transport. Note also that D{e) is discontinuous 
(but finite) at e = ±|A|. These discontinuities vanish 
for A = 0. On the other hand, if the Zeeman coupling 
dominates over the spin-orbit interaction, 



£R < |A| 



(7) 



the energy dispersion branches have stationary points 
only at zero wave vector. Here no van Hove singular- 
ity occurs in the density of states, apart from a step-like 
behaviour at minimum energy. These two cases er ^ || A| 
are illustrated in the left panel of Fig. [TJ Note that for 
er > |A| the avoided crossing at k — does not lead to a 



van Hove singularity since these stationary points occur 
at vanishing wave vector. 

In the case of an arbitrary magnetization direction, 
A = (A x , A y , A z ), a closer analysis shows that the sta- 
tionary points of the dispersion (fulfilling ds±/dk = 0) lie 
at wave vectors perpendicular to the in-plane projection 
of the magnetization, 



k- A = 0. 



(8) 



However, since this very general case does not seem to 
allow for further explicit results, we shall concentrate on 
a purely in-plane magnetization, A = (A x , A y ,0). Here 
the stationary points of the dispersion correspond to en- 
ergy minima and lie at wave vectors 



min 



± 



m*a (-A y ,A x) 



h 2 A 
The absolute minimum is given by 

£ rnin = £ -{Knin) = — ~^ E R ~ A 



(9) 



(10) 



(assuming a > 0), while another minimum occurs at k 



k tin with 



= e±(k 



+ 

min 
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A, 



(11) 



where the plus (minus) sign in e±(-) applies if A > er 
(A < er). These two cases are depicted in the right 
panel of Fig. [TJ The dispersion minima at finite wave 
vector k = fcjL„ are physically easily understood: In the 
case £ = e~ m , the Zeeman field the and spin coupling 
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provided by the Rashba interaction are parallel to each 
other leading to an energetic minimum for the appropri- 
ate spin direction, while for the higher minimum e = e^ nin 
these couplings are antiparallel. Note also that these en- 
ergetic minima remain at finite wave vector for arbitrar- 
ily large magnetization. Therefore, differently from the 
case of perpendicular magnetization, the van Hove sin- 
gularities in the density of states do not vanish for large 
Zeeman coupling. 



B. Discrete system 

The standard discretized version of the Hamiltonian 
([1]) on a square lattice reads 



— t y ] [ C m,n;a C "i+l,n:c 



V v = V /2 



V D = Vb/2 



C m.n:a C m,n+l;a + h.C.] 



m,n;a,/3 

+ ic t n,n ;a cr a, ; 3 C m+l,n;/3 + h.C] 

+ X! C ln,n;a^m.n ' V a ,/3C m ,n;l3 ■ (12) 
m,n;a,/3 

Here m and n label lattice locations with respect to the 
x- and y-axis, respectively, and a,/3 are spin indices. The 
hopping parameter t is related to the effective mass m* 
and the lattice constant a via t = ti 2 /2m*a 2 , and the 
parameter A is given by A = a /2a. To give a specific 
example, for a host material like gallium arsenide we have 
an effective band mass of m* = 0.06777io (with mo being 
the free electron mass) and a lattice spacing of a = h.QA 
leading to a hopping parameter t = 1.8eV. Typical values 
for the Rashba parameter are of order O.leVA such that 
we have typically A ~ O.Oli. 

For a spatially homogeneous impurity polarization of 
arbitrary direction, A mj „ = (Ax, Ay, A z ) the energy dis- 
persions are given by 

£±{k) = — 2t(cos(k x a) + cos(k y a)) 



± 



(A x - 2Asin(fc y a)) 5 



-{Ay + 2Xsm(k x a)) 2 + A 2 Z 



1/2 



(13) 



In order to make contact to the continuum model, one has 
to evaluate this dispersion for small wave vector, fco< 1, 
corresponding to the lower band edge, where it repro- 
duces the Eq. ^ up to a rigid shift of (— At) which is 
just half of the band width in the absence of magnetiza- 
tion. 

However, these dispersion relations (fT3"|) of the discrete 
system lead to very intricate conditions for stationary 
points which do not seem to be explicitly solvable. We 



FIG. 2: Four-terminal Hall bridge. The central region is 
described by the discretized Hamiltonian (|12|l incorporating 
Rashba spin-orbit coupling and magnetic impurities. 



therefore concentrate on the case of a homogeneous po- 
larization perpendicular to the plane, A m .„ = (0,0, A). 
Here the lower branch leads again to a singularity in the 
density of states provided that 

2A 2 >£|A|, (14) 

which is exactly the same as the condition ([3]). 



C. Hall bridge and Landauer-Biittiker formalism 

The above discretized system described by the Hamil- 
tonian (fl"2"l) is studied as the central region of a four- 
terminal Hall bridge shown in Fig. [5] In the follow- 
ing, we will investigate this system consisting of its cen- 
tral part and ideal semi-infinite leads without spin-orbit 
coupling and magnetization using the Landauer-Biittiker 
formalis m 25 ' 26 . We now briefly summarize the most im- 
portant features of this method as applied to the calcu- 
lation of the Hall conductivity. 

Within the Landauer-Biittiker approach, the Hall con- 
ductivity is given by 



h 



1 



1 



aH = V^~Vn = " 2 GuL + 2 GuR ' (15) 
where Ijj is the current flowing into lead up, and Gjjl and 
Gur are the conductances between lead up and left and 
between lead up and right, respectively. These quantities 
can be calculated via the following equation 



I P = Y,G pq (V p -V q ), 



(16) 



which describes the current flowing in lead p, where 

,2 



G v 



- T pq{ £ f) 



(17) 



4 



is the charge conductance between leads p and q. This 
quantity is proportional to the transmission function T pq 
defined as 



T pq = tv[T p G r T q G a ] 



(18) 



where the retarded (advanced) Green's function G r ^ 
and the line-width function V enter. These functions are 
the heart of the Landauer-Buttiker approach. Indeed, 
both Grecns's functions 



G r = [G a ] 



H 



and 



T q = i[E q - E«] 



(19) 



(20) 



depend on the retarded (advanced) self energy S r ( a ) . The 
retarded Y> r n „ and the advanced „ = T, r „\. self-energy 
terms describe the coupling between the two-dimensional 
electron gas (2DEG) in the central region and the four 
ideal scmi-infinitc leads and can be formulated as 



(21) 



where q r n ,, is the Green's function of the isolated semi- 
infinite lead. The only terms which enter into the self- 
energy are the coupling terms between each lead and the 
2DEG central region. The coupling between a lead and 
the central system give rise to a coupling matrix which is 
non-zero only for adjacent points i (lying at the edge of 
the 2DEG) and pi (lying at the leads edge which faces the 
2DEG). Since the lead Green's function can be evaluated 
analyticall y 25 ! 26 , the self-energy method can be exploited 
to deal with an infinite system, such as a semi-infinite 
lead, by calculating only the Green's function of a finite 
region. The analytical expression of our Green's function 
is derived in appendix [A"l 

Another important quantity which can be calculated 
within the Landauer-Buttiker formalism is the density of 
the states (DOS) given by 



D(e) 



1 



2-nN l a 



■M[A{e)] 



1 



2irN 2 w 



■tr[Im(G'')], (22) 



where A(e) = i[G r — G a ] is the spectral function, and 
N 2 is the number of lattice sites in the central region 
taken to be a square. The above expression is identical 
to another standard textbook result, 



D(s) 



1 



(27T) 2 
1 



d 2 k6(e-e(k)) 



dk- 



e(k)—e 



(23) 
(24) 



which can be used to derive , e.g., Eq.([S]). For an infinite 
2DEG with dispersion e(k) = —2t(cos(k x a) + cos(k y a)) 
the DOS can be calculated analytically leading to a loga- 
rithmic divergence at e — and a saturation to constant 



values at the edges of the band, i.e. D(e = ±4t) = j^ta 2 . 
On the other hand, in the case of a finite 2DEG central 
region and no coupling to the leads, the DOS is just given 
by a sum of S— peaks for values of Sf which match the 
eigenvalues of the Hamiltonian of the central conductor. 
As we shall see below, the anomalous Hall conductance 
is closely related to the DOS. 



III. RESULTS 

Let us now describe our numerical results based on the 
Landauer-Buttiker formalism outlined before. We first 
concentrate on disorder-free ballistic systems. 



A. Ballistic Hall Transport 

Here we present our results for a disorder-free cen- 
tral region with homogeneous magnetization A. We will 
both cover the case of magnetization perpendicular to 
the plane of the 2DEG, A = (0, 0, A) and the case of in- 
plane magnetization of various directions. In the latter 
scenario, the magnetization A can also be interpreted 
as a proper magnetic field B which, in a strictly two- 
dimensional system, couples only to the spin but not to 
the orbital degrees of freedom of charge carriers. 



1. Magnetization perpendicular to the 2DEG 

Let us first turn to the case of systems magnetized 
perpendicularly to the plane of the 2DEG. We consider 
a magnetization A = (0,0, A) and have evaluated the 
Hall conductance and the density of states according to 
Eqs. (fT5)) and (|22[) . respectively. As is must be, the Hall 
conductance vanishes for zero magnetization since the 
conductances Gul and Gur are identical and cancel out. 
With a finite Zecman coupling, however, a charge current 
Ijj starts to flow in lead UP signalling a finite Hall con- 
ductance. Fig[3] shows the Hall conductance along with 
the DOS for a Zeeman coupling A = 0.00U and a Rashba 
parameter of A = 0.0U as a function of Fermi energy 
Sf G [—At — A, 4t + A]. Both quantities plotted are per- 
fectly symmetric with respect to the band center Ef = 0. 
The Hall conductance is characterized by an oscillatory 
behavior over the entire energy range with particularly 
dominating peaks near the band edges; a smaller peak 
occurs also at the band center. On the other hand, the 
DOS shows the predicted logarithmically divergence at 
Sf =0 and in addition even more pronounced singulari- 
ties near the band edges at exactly the same positions as 
the peaks of the Hall conductance. 

As we wish to make contact to previous 
theoretical investigations on anomalous Hall 
transport in 2DEGs described by continuum 
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FIG. 3: (Color online) Hall conductance and DOS for a Zee- 
man coupling A = O.OOlt and a Rashba parameterA = O.Olt 
as a function for Fermi energy £/ £ [—At — A, At + A] and a 
linear system size of N = 30. The inset shows the behavior 
of the Hall conductance and the DOS close to the lower band 
edge. Both quantities are characterized by a simultaneous 
peak. 



model8 ll i l2^3 i l4 1 l5,l6 l l7 1 l8 l l9 1 20 l 2l 1 22,23 ) we wiU concen _ 

trate in the following on the peaks of Hall conductance 
and DOS close to the lower band edge. In FigfJ]we have 
plotted both quantities near the lower band edge for the 
same Rashba parameter and system size as before, but 
different Zeeman couplings A. As seen in the figure, 
extrema of the DOS and the Hall conductance occur 
at the same position in energy, independently of the 
regime of Zeeman coupling. The close correspondence 
between extrema of the DOS and the Hall conductance 
will be an important finding for our further analysis of 
anomalous Hall transport. Note that the maxima of the 
DOS at lower energies become weaker with increasing 
Zeeman coupling A. Such a behavior can indeed be 
expected from the analytically accessible properties of 
the infinite system discussed in section [TTJ Here the 
divergent van Hove singularities in the DOS vanish if 
the magnetization dominates the spin-orbit coupling. 
The systems investigated in this work numerically 
are obviously different as they are finite and coupled 
to semi-infinite leads lacking spin-orbit interaction. 
However, the above observations shall still guide our 
intuition regarding the interplay between magnetization 
and spin-orbit coupling. 

Fig. [5] displays the Hall conductance near the bottom 
of the band for a Rashba parameter of A = O.Qlt and 
various Zeeman couplings A with the linear system size 
varying from N = 30 to N = 50. For small magneti- 
zation A < O.OOlt the height of the Hall conductance 
peaks is approximately independent of the system size, 
while for larger Zeeman couplings 0.001 < A < 0.0075i 
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FIG. 4: (Color online) DOS (red dashed line) in units 1/ta 2 
and Hall conductance (blue solid line) in units of e 2 /h near 
the lower band edge for the same system size and Rashba 
coupling as in Fig. [3] but various Zeeman couplings A. In all 
cases an obvious correspondence between the extrema DOS 
and Hall conductance occurs. 



slightly grows with increasing system size. For even 
larger A > 0.0075i a decrease is observed for large sys- 
tem sizes. This qualitative behavior persists in a range 
of Rashba parameters A G [0.005i, 0.02i] with the above 
threshold values for the Zeeman coupling A being ap- 
proximately unchanged. While our above finite-size data 
for the height of the Hall conductance peaks does not 
seem to allow for a unambiguous extrapolation to the 
thermodynamic limit, the suppression of the Hall con- 
ductance at large Zeeman splittings A > 0.0075i is cer- 
tainly consistent with the analytical observations in the 
infinite system. At large Zeeman couplings, the singular- 
ity in the DOS close to the band edge disappears, and, 
in turn, the Hall conductance vanishes. 

Moreover, as also seen in Fig. [51 for all Zeeman cou- 
plings, the position e* N of the peak shifts to lower en- 
ergies , i.e. towards the bottom of the band, with in- 
creasing system size. Fig. [5] shows the finite-size be- 
havior of the peak position for same six data sets as in 
Fig. All data sets can be smoothly fitted by an ex- 
ponential function which allows for an extrapolation to 
the limit of an infinite system, e* = limjv^oo 

The 

dependence of e* on the magnetization A is shown in 
Fig. [71 Clearly, e* linearly decreases with increasing A. 
Note, however, that for large Zeeman splitting the Hall 
conductance peaks arc suppressed with increasing sys- 
tem size although their position can still be followed as a 
function of N. On the other hand, even for the smallest 
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FIG. 5: (Color online) Hall Conductance near the bottom of 
the band for a Rashba parameter of A = 0.0 It and various 
Zeeman couplings A. The linear system size varies from N = 
30 (blue, right), N = 35 (red), TV = 40 (green), N = 45 
(yellow), and N = 50 (black, left). 



Zeeman gap of A = O.OOli considered here, the infinite- 
volume peak position e* lies at an energy where only the 
lower dispersion branch of Eq. (fl3|) is occupied. This is in 
accordance with recent theoretical predictions based on 
diagrammatic perturbation theory by Nunner et al. who 
concluded that a finite Hall conductance can only occur 
at low Fermi energies such that only the lower subband 
is occupie d 18 ' 20 . 



2. In-plane Magnetization 
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FIG. 6: (Color online) Dependence of the position e* N of the 
Hall conductance peak on the system size for same data sets 
as in Fig. [5] The Rashba parameter is A = O.Olf. 
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FIG. 7: Energy position e* of the Hall conductance peak at 
the lower band edge extrapolated to an infinite system versus 
the magnetic coupling A: the dependence is well fitted by a 
linear function. 



We now turn to the case of magnetic impurities polar- 
ized in the plane of the 2DEG. As shown in Fig. [21 the 
Hall current is measured along the y-direction in leads 
up and down, while the Hall Voltage is applied along 
the x-dircction between leads left and right. Here we do 
not observe a Hall current if the in-plane magnetization 
is parallel to the voltage since here for charge carriers 
with wave vector along the y-direction the spin coupling 
resulting from the Rashba interaction and the Zeeman 
coupling are just parallel to each other. However, a finite 
Hall current can occur for other in-plane directions of the 
magnetization. 

Fig. [5] shows the Hall conductance for several direc- 
tions of magnetization. For magnetization along the 
x-direction no Hall current occurs, and for a magneti- 



zation pointing in the y-direction with A y = 0.0075i 
we observe a maximum of the Hall conductance which 
occurs at the same energy as for the previous case of 
strictly out-of-plane magnetization A 2 = 0.0075t, but 
is smaller in magnitude. For a magnetization pointing 
in the (1, 1, 0)-direction with A x = A y = 0.0053t (ful- 
filling |A| — 0.0075i) we find a conductance peak at 
the same position in energy but of different shape. If 
the magnitude of the Zeeman splitting is increased to 
A x = A y = 0.0075£ the peak approximately maintains 
its shape but is shifted towards lower energies. We fo- 
cus now on an in-plane magnetization totally polarized 
along the y-direction. Fig. [9] in the "in-plnae analogue" 
of Fig. [Jf] and shows the correspondence between the ex- 
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FIG. 8: (Color online) Hall conductance for a Rashba pa- 
rameter of A = O.Olt and various directions of magnetization 
with the non-zero components of A stated in the legend. For 
magnetization along the x-direction no Hall current occurs. 
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FIG. 9: (Color online) DOS (red dashed line) in units 1/ta 2 
and Hall conductance (blue solid line) in units of e 2 /h near 
the lower band edge for an in-plane magnetization polarized 
along the y-direction, for a system of linear size N = 30 and 
for three different magnitudes of the Zeeman splitting. This 
picture is analogous to Fig. U for the case of a perpendicular 
magnetization. Even in the case of an in-plane magnetization 
we observe the correspondence between the extrema DOS and 
Hall conductance. 
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FIG. 10: (Color online) Hall Conductance near the bottom of 
the band for a Rashba parameter of A = O.Olt in the case of an 
in-plane magnetization along the y-direction and for various 
Zeeman couplings A y . The linear system size varies from 
N — 30 (blue, right), N = 35 (red), N = 40 (green), N = 
45 (yellow), and iV = 50 (black, left). This picture can be 
directly compared with Fig. [5] for the case of a perpendicular 
magnetization. 



trema of the DOS and those of the Hall conductance for 
an in-plane magnetization totally magnetized along the 
y-direction for a system of linear size N = 30 and for 
three different values of the magnitude of the Zeeman 
splitting. As in the case of a perpendicular magnetiza- 
tion, see Fig. [4j we observe a correspondence between 
the extrema of the DOS and the Hall conductance at the 
same Fermi energy. 

In Fig. [TU] we show the Hall conductance for an in- 
plane magnetization along the y-dircction varying be- 
tween A y = O.OOlt and A y = O.Olt, and linear system 
sizes between N — 30 and N = 50. The Rashba pa- 
rameter is again A = O.Olt. Obviously, the position of 
Hall peaks shifts to lower energy with increasing systems 
size, analogously as in Fig. [5] for the case of a perpendic- 
ular magnetization. Moreover, for the smallest Zeeman 
coupling considered here, A y = O.OOlt, the height of the 
peaks clealy grows with systems size, while for the largest 
Zeeman coupling of A y = O.Olt the opposite behavior is 
observed. However, we cannot outline any trend for in- 
termediate values of the Zeeman coupling, see Fig. [TU] 
b), c), d), e). Another interesting finding is that the en- 
ergetic position of the Hall signal coincides for in-plane 
and perpendicular magnetization of the same magnitude. 
This result is shown in Fig. 111! where the linear size of 
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FIG. 11: (Color online) Hall conductance and DOS for a sys- 
tem of linear dimension N = 40 and for a perpendicular mag- 
netization of amplitude A z = 0.014 (blue dashed line) and an 
in-plane magnetization A y = O.Oli (green solid line) 



the system is N = 40 and the Zeeman coupling is cho- 
sen to be A z = A y = O.Olt. Here, in the two cases of 
perpendicular and in-plane magnetization the shape of 
the DOS and Hall conductance peaks varies but their 
maxima exactly coincide . The dependence of the Hall 
conductance on direction and magnitude of the in-plane 
magnetization is to be explored further in the future. 



B. Magnetically disordered systems 

So far we have studied homogeneously magnetized sys- 
tems with each lattice site carrying a magnetic ion whose 
spin provides a Zeeman field for the carrier spins. We now 
consider the case where only a given fraction v of lattice 
sites is occupied by a magnetic ion. This scenario ac- 
counts for the situation in real ferromagnetic semiconduc- 
tor nanostructures. We will concentrate again on magne- 
tizations along the z-direction perpendicular to the plane 
of the 2DEG. To be specific, we choose at random a given 
fraction of lattice sites to be occupied with a magnetic 
ion and average our results for the Hall conductivity over 
typically 20 of such disorder realizations, which, by in- 
spection of the data, turns out to be sufficient. In order 
to facilitate the comparison with our previous results for 
magnetically homogeneous systems we adjust the magni- 
tude of each local coupling A mi „ such that the average 
Zeeman coupling A := (J2 m n A m , n )/N 2 is constant, i.e. 
A mj „ = Ajv. Fig. [P2l shows the Hall conductance for 
different fractions of magnetically occupied sites for a 
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FIG. 12: (Color online) Hall conductance close to the lower 
band edge for different fractions of magnetically occupied sites 
for a Rashba parameter of A = 0.014 and an average magneti- 
zation of A = O.OOlt. The linear size of the system is N=40. 



Rashba parameter of A = O.Olt and an average magne- 
tization of A = O.OOlt. The data is averaged over 20 
random disorder configurations. Fig. [T3] shows the de- 
pendence of the Hall conductance both on the fraction 
of magnetized sites as well as on system size. Again, 
the energetic position of the hall maxima moves towards 
lower energies with increasing system size, cf. Figs. [5TT0l 
Moreover, for fractions v = 0.3 . . . 0.9 of magnetized sites, 
the position of the Hall peaks is very close to that of the 
uniform magnetization, v = 1, and the height of peaks 
is approximately constant in all cases. This is different 
from the smallest fraction considered here, v = 0.1, where 
height increases with system size, and also the energetic 
positions differ from the other cases. 



IV. CONCLUSIONS AND OUTLOOK 

Anomalous Hall transport is still the subject of a long- 
standing theoretical discussion. In the present paper we 
have report on the, to our knowledge, first investigation 
of this phenomenon using the Landauer-Biittiker formal- 
ism. Specifically, we have studied a two-dimensional elec- 
tron gas with Rashba spin-orbit coupling and a magne- 
tization provided by localized magnetic moments. Our 
system is described by a discretized tight-binding model 
in a four-terminal geometry. We have considered both 
the case of ballistic transport in disorder-free systems 
with homogeneous magnetization as well as magnetically 
disordered systems. In the former case we have also 
distinguished between different directions of the magne- 
tization. In particular, a magnetization lying entirely 
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FIG. 13: (Color online) Hall conductance maxima for systems 
of linear dimension N = 30 (blue, right), N = 40 (green), 
N — 50 (black, left) and for several fraction of magnetically 
occupied sites. 



APPENDIX A: GREENS FUNCTION OF A 
SEMI-INFINITE LEAD 

Here we give some more technical details regarding the 
calculation of the Greens function of Eq. (j2"Tj) for a semi- 
infinite lead. 

For a semi-infinite non-interacting lead with hard wall 
confinement and a constant width L, the transverse wave 
functions arc 



Xm(y) 



sin 



or 



Xm {Uj 



2 . ( mnj 
■ sin 



M + l \M + 1 



(Al) 



(A2) 



where yi = ja and M is the number of sites in the 
transverse direction, such that L = (M + l)a. The lon- 
gitudinal lattice wave functions are 



<t>k(x) 



sin(A;ar), 



(A3) 



which, by substituting x = a, which means that we are 
considering points x at the first slice of the semi-infinite 
lead, transforms into 



(f> k (x = a) 



■ sin(fca). 



(A4) 



in the plane of the 2DEG can also be interpreted as a 
genuine magnetic field which couples, in a strictly two- 
dimensional system, only to the spin of the electron but 
not to its orbital degrees of freedom. 

In particular, we have demonstrated numerically a 
close connection between singularities in the density of 
states and peaks in the Hall conductance close to the 
lower band edge. Consistent with previous theoretical 
studies based on diagrammatic perturbation expansions, 
these peaks occur at Fermi energies where only the lower 
dispersion branch is occupie d 18 i 20 . Moreover, for large 
magnetization the Hall conductance is, along with the 
density of states, suppressed. This numerical finding can 
be understood from analytical properties of the underly- 
ing model in the limit of an infinite system. 

Future investigation will include a more detailed un- 
derstanding of anomalous Hall transport in the presence 
of an in-plane magnetization, and the effects of magnetic 
disorder. 
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Finally, at the first slice x = a the total wave function 
reads 



4>m,k 



-Xm{y) sin(fca), 



(A5) 



which is an eigenfunction of the Hamiltonian -ffoV J m,fc = 
E m ,k4>m,k, where the dispersion relation is given by 



E,, 



2t 1 - cos 



M + l 



+ 2t(l - cos(fca)). (A6) 



Inserting the wave function in the eigenfunctions ex- 
pression of the Greens function, see&, the Greens func- 
tion at the first slice x = a is 



G r (a,y t ;a,y j ) 



2a 



E 

m,k>0 



Xm(yi)x* m (yj)sm 2 ( ka ) 

E - E m ,k + if] 



(A7) 

In the limit L — > oo, we may replace the sum over k 

by the integral (L/n) J? dk and substitute ka = 6 to 
obtain 

G r ( a ,y^a,yj) = l/irfEmXmiVijXmiVo) /Aa\ 
/„" d6 sin 2 e/(Q + cos + irj), 1 ' 



where we made the following replacement 

E mir 

Q = 2 + cos — . 

^ 2t M + l 



(A9) 



10 



In Eq. (|A8[) the integrand function sin 2 9/(Q+cos 9+irj) is 
an even function of 9, therefore the integral can be writ- 
ten as a symmetric integral from — 7r to it. By writing 
in the integral the sine and cosine functions in terms of 
their exponential form and performing the substitution 
z = exp(i9), the integral turns into a closed contour inte- 
gral along the unit circle in the complex plane, allowing 
us to write the Green function as 

G r (a,y l ;a,y j ) = l/(2mt) J2 m Xm{yi)Xm{yj) /» 1n v 
J c dz(l-z 2 )/(z 2 + 2zQ + l), [AW) 

where J c stands for the integration on a closed circuit of 
radius | z |= 1. By solving the integral with the use of 
the theorem of residues, we obtain 

G r (a,y l ;a,y J ) = l/(2mt) £ ro Xm{Vi)Xm{Vj) / A11 x 
2™R 2o (l-z 2 )/(z 2 + 2zQ + l), ^ Allj 

where R Zo indicates the residual calculated at the pole 



z = zq which depends on Q, i.e. 



20 



Q + V Q^l for Q > 1 
Q - y ^ 2 ^! for Q < -1 (A12) 
-Q + zv/I - Q2 for | Q |< i 



Inserting the pole in Eq. (|Allj ). wc obtain the final ex- 
pression for the green function 

G r (a, yi ;a, yj ) = - £ Xm{Vi)x* m {Vi)F{Q), (A13) 



with 



Q-VQ^Zl for Q>1 
F(Q) = { Q+ v/^FZL for Q < - 1 ( A14 ) 
Q - W l - Q 2 for I Q l< i 
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